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Motivation: Sign problem, Silver Blaze problem 


Yuya Tanizaki (RIKEN, RBRC) 


Lefschetz-thimble path integral 


Sign problem of finite-density QCD 


QCD 
Fundamental theory for quarks and gluons 


Neutron star 


e Cold and dense nuclear matter 


@ 2Msun neutron star (2010) 


Neutron star merger 
e Gravitational-wave observations (2016) (image from NASA) 


Reliable theoretical approach to equation of state must be developed! 


Zaco(T. 1) = [ DA Det(P(A, pa) +m) exp C Svu(). 


quark gluon 


Sign problem: Det(D(A, uq) 4- m) Z 0 at uq z 0. 
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Sign problem of finite-density QCD 
QCD & |QCD| 


ZQCD = [ra (det D) e SYM | Ziacp| = [ea [det P| eM, 


If these two were sufficiently similar, we have no practical problems. 
However, it was observed in lattice QCD simulation that at T = 0 
(e.g., Barbour et. al. (PRD 56 (1998) 7063)) 


NB QCD 


E > 
-mr/2 -myw/3 Ha 
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Baryon Silver Blaze problem 


The curious incident of the dog in the night-time (Holmes, Silver Blaze). 
Problem: Show that ng = 0 for ui, < my/3 using path integral. 
(Cohen, PRL 91 (2003) 222001) 


Current situation: For jug < m; /2, the problem is solved. 
Quark det. is 


Det(P(A, jg) +m) P" 
Det(D(A, 0) +m) Il Greve), 


0O<Re(AA)<fkg 


and ess-min4(Re(A4)) = m, /2. (Cohen, PRL 91 (2003) 222001, 
Adams, PRD 70 (2004) 045002, Nagata et. al., PTEP 2012 01A103). 
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Spectrum of ^4(I) 4 + m) 


Gap of {An}n = Spec (^4(D 4 + m)} gives the pion mass (Gibbs). 
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Lattice study of the quark spectrum and the Dirac determinant. 
(Nagata et. al., PTEP 2012 01A103) 
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Method: Path integral on Lefschetz thimbles 
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Lefschetz-thimble path integral 


Sign problem of path integrals 


Consider the path integral: 


ja ] «eec sta 


e S[x] is real = No sign problem. Monte Carlo works. 
e S|[x] is complex = Sign problem appears! 


If S[x] € C, eom S'[x| = 0 may have no real solutions z(t) € R. 


Idea: Complexify x(t) € C! 
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Lefschetz thimble for Airy integral 
Airy integral is given as 


3 
Ai(a) = [ri (= + 2 


Complexify the integration variable: z = x + iy. 


4 


Integrand on R, and on A 
(a — 1) 
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Rewrite the Airy integral 
There exists two Lefschetz thimbles J, (c = 1,2) for the Airy 


integral: 
Ai(a) Son / “an am 
= ya — — az]. 
= J, 2T p 3 


ng: intersection number of the steepest ascent contour Ke and R. 
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Lefschetz decomposition formula 


Oscillatory integrals with many variables can be evaluated using the 
“steepest descent” cycles 7,: (classical eom S’(z,) = 0) 


/ d”gz e- 560 = 3s R) | d”z e80, 


o 


Jo are called Lefschetz thimbles, and Im[S] is constant on it: 


ma =}, dzi(t) _ eS 


t>—00 dt Oz? 


I= 20) 


(Ko, R): intersection numbers of duals K, and R” 
(Ke = {2(0)|2(c0) = #2}. 

[Witten, arXiv:1001.2933, 1009.6032] 

[Christoforetti et al. (PRD(2012)), Fujii et al. (JHEP 1310), etc.] 
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Analysis: Semi-classical analysis of the one-site Hubbard model 
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Lefschetz-thimble path integral 


One-site Fermi Hubbard model 
One-site Hubbard model: 


A = Uff, — u(fu + f). 


Fock state gives the number density immediately: 


wA 


(YT, Hidaka, Hayata, 1509.07146)(cf. Monte Carlo with 1-thimble approx. gives a wrong result: 


Fujii, Kamata, Kikukawa,1509.08176, 1509.09141; Alexandru, Basar, Bedaque,1510.03258.) 
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Path integral for one-site model 
Effective Lagrangian of the one-site Hubbard model: 


| P(r) 
m 2U 


"dr — (U/2 + i¢(T) + u)] v. 


The path-integral expression is (p = feo r)dr/ B) 


B / : 2. i 
- du (1 Bligt+p+U/2) By" /2U | 
Von re 


Fermion Det 


Integrand has complex phases causing the sign problem. 


ọ is an auxiliary field for the number density: 
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Sign problem and fermion determinant 


One-site Hubbard model: 
Uae —B(—U/2-1) oib)? 
Det |0, — jm enit = (1 Fe MER), 


Quark determinant in QCD: 


Det UP, +m) — u] = N(A) I] (14-678 6;7u-i6)) (]4e-Blertutien)), 


£j20 
where the spectrum of ^4(ID , + m) is 
Atıın) = €j (A) = id; (A) + (2n + Du. 


Minimal value of £(A) = m; /2. 
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Silver Blaze problem for u < —U/2, u < m4/2 


One-site Hubbard model: As BU >> 1 and —U/2 — u > 0, 


U 
Det lo = (» + 7 + e) = (1 + EURER)" ~]. 


The sign problem almost disappears, so that J, ~ R. 
Finite-density QCD: As 6 — oo and u < m,/2, 


Det ud, +m) — ul 
Det l4 D 4 F m)| 


The sign problem disappears by the reweighting method. 
=> Lefschetz thimbles ~ Original integration regions 


> 1 
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Analysis 


Flows at u/U < —0.5 (and u/U > 1/5) 
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Figure: Flow at u/U=-1. Ix 


-0.25 


-0.5 


cR. 


; dz e- 909, 


Number density: n, = 0 for u/U < —0.5, n, = 2 for u/U > 1.5. 


(YT, Hidaka, Hayata, 1509.07146) 


Z- 
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Silver Blaze problem for u > —U/2, u > m,/2 


One-site Hubbard model: At each real config., the magnitude is 
exponentially large: 


Det lo — (» + 5 + e) SO (AUWA 


This large contributions must be canceled exactly in order for n = 0. 


Finite density QCD: The situation is almost the same, since 


Det(D(A, uj) tm) _ 
Det(D(A,0)+m) — 


I 9»56-24. 


Re(AA) «uq 


but ng = 0 for u, S my/3. 
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Analysis 


Flows at —0.5 < u/U < 1.5 
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Figure: Flow at u/U = 0 


-0.5 


-0.75 


py p 1. 


This value is far away from n = Im (z)/U=0, 1, or 2. 


m 


Complex saddle points lie on Im(z,,)/U 


20 / 33 


Lefschetz-thimble path integral Oct 6, 2016 0 SBU 


Yuya Tanizaki (RIKEN, RBRC) 


Curious incident of n in one-site Hubbard model 


We have a big difference bet. the exact result and naive expectation: 


95 0 0.5 1 15 ajo 


This is similar to what happens for QCD and |QCD|. 
p/U = -0.5 & py =m, /2. 
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Complex classical solutions 


If GU > 1, the classical sol. 
for —0.5 < u/U < 1.5 
are labeled by m € Z: 


u[U x 1/2 


U 
Ea cR By +2rmT. 


At these solutions, the classical actions become 


Sox -£ (+, 
Re (Sin — So) ~ 2m m?, 
Im S. c 2mm (5 4- 1). 
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Semiclassical partition function 


Using complex classical solutions Zm, let us calculate 


= 1 
Za ‘= p» em 2 e 50009, (s (5 J 5) uud . 


This expression is valid for —1/2 S u/U S 3/2. 
YANI 


-1/2 0 1 3/2 
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Important interference among multiple thimbles 


Let us consider a “phase-quenched” multi-thimble approximation: 


Za) =) le | = AA 


e Lee-Yang zeros cannot appear at u/U = 0, 1. 
e One-thimble, or “phase-quenched”, result: n ~ u/U + 1/2. 


Consequence 


To understand the Silver Blaze problem, we need interference of 
complex phases among different Lefschetz thimbles. 


(cf. Particle Productions: Dumulu, Dunne, PRL 104 250402) 
(cf. Hidden Topological Angles: Behtash, Sulejmanpasic, Schafer, Unsal, PRL 115 041601) 
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Numerical results 


Results for BU = 30: (1, 3, 5-thimble approx.: Jo, Jo U Ji, and Jo U J41 U Ja ) 
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Necessary number of Lefschetz thimbles ~ BU/ (27). 


(YT, Hidaka, Hayata, NJP 18 (2016) 033002, arXiv:1509.07146[hep-th]) 
Yuya Tanizaki (RIKEN, RBRC) 
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Bonus 


Bonus: Complex Langevin method 
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Bonus 


Complex Langevin method 


Complex Langevin has been regarded as a sign-problem solver via 
stochastic quantization (Klauder, PRA 29, 2036 (1984), Parisi, PLB 131, 393 (1983)): 


dz,(0) 0S 
= — 5 00) + VANO). 


0: Stochastic time, n: Random force (n(0)n(0')),, = 20(0 — 0"). 


Properties: 


e Numerical cost is very cheap. 

e Ito calculus shows (O(z,(oo))); solve the Dyson-Schwinger eq. 
e Sign problem does not appear. 

e But, it fails in some cases. = When does it fail? 
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Bonus 


Semiclassical incorrectness of CL method 


If A is small enough, we can show a sufficient condition for incorrect 


behaviors of CL method. 
(Hayata, Hidaka, YT, 1511.02437) 


Since A < 1, CL distribution would accumulate around {zo}: 


"Op 0) wt oO es S esO(as). 
Assume for contradiction that CL method is correct, then 
(Oley + f dee SM0(@) 
R 


1 
z 2C. R") / dz e-5/^9(;), 
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Bonus 


Semiclassical inconsistency 


In the semiclassical analysis, one now obtains (for dominant saddle 
points) 


- (Ko, R”) Ith e $69)/h. 
Z S" (zo) 


The right hand side can be complex, which contradicts with c, > 0! 
(Hayata, Hidaka, YT, 1511.02437) 


We show that the complex Langevin is wrong if 
e There exist several dominantly contributing saddle points, and 
e Those saddle points have different complex phases. 
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Bonus 


Proposal for modification 
Assume as a working hypothesis that 


_ (Ko, R) 2mh SG 


S"(z,) 


Because of the localization of probability distribution P, it would be 
given as 


P- ` CoP,, supp(P,)risupp(P,) = 0. 


o 


Assumption means “CL = phase quenched multi-thimble approx.” : 


(Ko, IR) 


2rh 
S" Za) 


eh 


Oz. ). 


(Olha = 


[02 
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Bonus 


Proposal for modification (conti.) 


If so, defining the phase function 


|S" (ze)| —iIm S(z;,)/h = 
S" (Ze) ii id Xsupp(Ps) (2, Z), 


we can compute 


(©(2,, Zn)O(Zn))n 
((25, 24))q l 


This new one is now consistent within the semiclassical analysis. 
(Hayata, Hidaka, YT, arXiv:1511.02437[hep-lat]) 

Caution: Our proposal evades inconsistency, but is not necessarily 
correct. Can we improve the proposal? 


Cej 
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Bonus 


Complex Langevin study of one-site fermion model 


BU =30 p/U=0 25 
= --- exact 
20l H= w/o phase 

= K w/ phase 
S 15 
+ T 
E = 10 -— 
| 
KD 0.5 
E 
«Y 
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(Hayata, Hidaka, YT, 1511.02437) 


Complex Langevin method cannot study the Silver Blaze problem of 
this model without introducing the reweighting procedure. 
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Summary & Conclusion 


Summary and Conclusion 


e Picard-Lefschetz theory gives a suitable framework for 
saddle-point analysis even if S(@) takes complex values. 


e One-site Hubbard model is a nice toy model to play with the 
sign problem. 


e Destructive and constructive interference of complex phases 
among Lefschetz thimbles play a pivotal role for the baryon 
Silver Blaze. 
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